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Much of the present day qualitative phenomenology of the fractional quantum Hall effect can be 
understood by neglecting the interactions between composite fermions altogether. For example the 
fractional quantum Hall effect at v = n/(2pn ± 1) corresponds to filled composite-fermion Landau 
levels, and the compressible state at v — l/2p to the Fermi sea of composite fermions. Away from 
these filling factors, the residual interactions between composite fermions will determine the nature 
of the ground state. In this article, a model is constructed for the residual interaction between 
composite fermions, and various possible states are considered in a variational approach. Our study 
suggests formation of composite-fermion stripes, bubble crystals, as well as fractional quantum Hall 
states for appropriate situations. 
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I. INTRODUCTION 

In two dimensions, electrons subjected to a strong 
magnetic field avoid one another by capturing an even 
number of quantized vortices and turning into compos- 
ite fermions [1,2]. The composite fermions interact much 
more weakly than electrons. In fact, tremendous progress 
can be made towards understanding the dramatic be- 
havior of this system by treating composite fermions as 
non-interacting. It is not that there is no interaction be- 
tween the composite fermions, but it is weak and often 
does not alter the nature of the state in a qualitative 
manner, i.e., can be treated perturbatively. In such sit- 
uations, while the inter-composite fermion (CF) interac- 
tion is surely important for quantitative considerations, 
it does not affect the qualitative phenomenology, and can 
be neglected altogether when the aim is to describe the 
qualitative phenomenology. 

For example, the odd-denominator fractions v = 
n/(2pn ± 1) map into v* — n of composite fermions car- 
rying 2p vortices. Here, the system has a gap even in the 
absence of inter-CF interactions - in stark contrast to the 
electron system at v = n/(2pn ± 1) which had an enor- 
mous degeneracy in the absence of interactions - and it 
is possible that for many of these filling factors, the gap 
would not disappear as the inter-CF interaction is slowly 
turned on to its physical value. If that were the case, 
the inter-CF interaction is unimportant at a qualitative 
level. This explains the fractional quantum Hall effect [3] 
(FQHE) at filling factors v = nj (2pn± 1) in terms of the 
integral QHE [4] of composite fermions [1]. 

Another example is at fractions v = l/2p, which map 
into composite fermions at zero effective magnetic field. 
If the interactions between composite fermion are ne- 
glected, a Fermi sea is obtained [5]. Conceptually, it is 
more difficult to justify the neglect of inter-CF interac- 



tions in this there is no gap for the Fermi sea. 

Nonetheless, at least there is a well defined ground state 
for composite fermions, unlike for non-interacting elec- 
trons at v = l/2p, and various experimental studies [6] 
have shown that the CF Fermi sea is a good starting 
point for many purposes. 

We recall that at relatively high temperatures and low 
mobilities only the integral quantum Hall effect was seen, 
which is explained by neglecting the interaction between 
electrons. As the sample quality improved and tempera- 
ture lowered, FQHE was observed, which is understood 
in terms of weakly interacting composite fermions. While 
the model in which the composite fermions are taken as 
non-interacting has been strikingly successful, one may 
ask if another set of new structures would emerge as the 
experiments further improve. This question has moti- 
vated our investigation of the subtle physics arising from 
inter-CF interactions. Fortunately, there already exists 
a good microscopic approach for investigating this issue, 
in the form of the wave functions for composite fermions 
[1]. Even though these wave functions are motivated by 
the physics of independent composite fermions, they give 
an excellent description of the interaction between the 
composite fermions [7,8,2]. In fact, the wave functions 
practically give the exact solution to the problem, which 
implies that they incorporate the full interaction effects. 
For example, consider the low-energy branch of excita- 
tions, described as an exciton of composite fermions. If 
the composite fermions were non-interacting, this branch 
would not have any dispersion. In reality, there are oscil- 
lations in the energy as a function of the wave vector, 
arising from the residual interaction between the CF- 
particle and the CF-hole. The dispersion computed from 
the wave functions provides an accurate account of the 
oscillations. The wave functions similarly obtain the en- 
ergy of two CF particles or two CF holes quite accurately. 
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The inter-CF interactions are of course always relevant 
to quantitative issues. However, there are two situations 
in which they can also make a qualitative difference in 
the physics. The first situation is when the inter-CF in- 
teraction overwhelms the gap of the non-interacting CF 
model, thereby destroying the FQHE [9,10]. It has been 
shown that for v < 1/9 as well as for most fractions 
in higher Landau levels, the FQHE undergoes an exci- 
tonic instability, i.e., the energy of the exciton becomes 
smaller than the energy of the uniform filled CF-Landau 
level (LL) state. The CF theory thus not only tells us 
where FQHE occurs but also where it does not. (Recent 
experiments [11] have shown evidence for FQHE at 1/9 
at finite temperatures. Even though this observation is 
not presently understood, it ought to be stressed that it 
is not inconsistent with the earlier theoretical predictions 
which only address the zero temperature phase.) A pair- 
ing instability has also been suggested at v = 5/2, i.e., at 
the half filled second LL, where the residual interaction 
between the composite fermions is attractive [12,13]. The 
resulting instability of the CF Fermi sea causes a gap to 
open up here, producing a plausible scenario for the 5/2 
fractional quantum Hall effect [14]. 

The other situation where the inter-CF interaction is 
expected to be important is away from the special filling 
fractions v = n/(2pn± 1), where the topmost CF-LL is 
only partially occupied, and the ground state has a large 
degeneracy if the inter-CF interaction were to be turned 
off. The degeneracy for independent composite fermions 
is much smaller than for independent electrons even here, 
making the former a better defined starting point. The 
weak residual interaction between the composite fermions 
will decide the nature of the true ground state here. The 
present article deals with such situations. 

The first possibility that comes to mind when compos- 
ite fermions partially occupy a Landau level is that of 
their fractional QHE. This could happen when the com- 
posite fermion filling factor is v* = n ± n/(2ph ± 1), 
which correspond, as usual, to electrons filling factors 
given by v = v* /(2pv* ± 1). We will consider the spe- 
cial case when the base composite fermions are carrying 
two vortices, denoted by 2 CFs, and the 2 CFs in the top- 
most partially occupied CF-LL capture 2 additional vor- 
tices to become CFs (composite fermions carrying four 
vortices) and condense into an incompressible state with 
n filled Landau levels. The full state is then made up 
of two flavors of composite fermions, with 2 CFs in the 
lower filled CF LLs, and 4 CFs in the topmost CF-LL. 
Such states are quite analogous to the FQHE state at 
v = 2 + 1/3, where electrons in the second LL convert 
into composite fermions whereas the electrons in the low- 
est LL remain electrons. Such states have been consid- 
ered in the past [1,15,16], dubbed mixed flavor states of 
composite fermions. Of course, depending on the filling, 
it is also possible that composite fermions would capture 
additional vortices to form a Fermi sea or a paired state. 



Yet another possibility is for them to not capture any 
new vortices but rather to phase separate in some man- 
ner, possibly through formation of stripes. Our aim in 
this article is to consider these possibilities as quantita- 
tively as presently possible over a range of filling factors. 

It should be emphasized that while these states are 
easily viewed in terms of composite fermions, they are 
fantastically complicated in terms of electrons. Consider 
the possibility of pairing at v = 3/8, for example [16]. 
First, all electrons capture vortices to turn into 2 CFs at 
v* = 1 + 1/2; then, those in the half filled CF-LL put 
on two more vortices, and attempt to make a 4 CF Fermi 
sea; but the Fermi sea is unstable to pairing; pairing 
opens up a gap to produce fractional quantum Hall effect. 
Such a state could perhaps not be envisioned without the 
knowledge of composite fermions. 

In an earlier short article [17], we had studied compos- 
ite fermion fillings given by v* = n + 1/2 by this method. 
The present article gives further details and also consid- 
ers a broader range of filling factors. The plan of the 
article is as follows. In the next section, we enumerate 
the various approximations that go into the calculation. 
Sect. Ill discusses how we obtain the pseudopotentials 
for the inter-CF interaction, and how we obtain from it a 
real space interaction between composite fermions. Cer- 
tain preliminary observations are made in Section IV. 
The variational wave functions are given in Section V, 
and Section VI discusses the results of our calculations. 
In Section VII, we briefly investigate the effect of finite 
thickness on our results. 

II. ASSUMPTIONS 

We will be studying composite fermions at fillings of 
the type v* = n + is*, where n is the number of filled com- 
posite fermion Landau levels and v* is the filling fraction 
of the topmost composite fermion Landau level. In order 
to make progress, we will make several assumptions. 

The system will be taken to be fully spin polarized, 
which is a valid approximation for sufficiently high mag- 
netic fields. The method can be adapted to situations for 
which the states are only partially spin polarized, [16] but 
here we will confine our attention to fully spin polarized 
systems. 

We will map below the problem of interacting elec- 
trons into that of interacting composite fermions only of 
the topmost partially filled CF Landau level. The in- 
teraction between composite fermions will be obtained 
by "integrating out" the composite fermions of the lower 
filled CF-LLs. 

We will assume that the filled CF Landau levels can be 
treated as inert and the problem can be defined solely in 
terms of the composite fermions in the topmost partially 
filled CF Landau levels. This should be qualitatively 
correct between well defined FQHE states described by 
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filled CF Landau levels. It should also be a quantitatively 
accurate approximation when the CF-LL separation is 
large compared to the inter-CF interaction. Indeed, for 
filling factor in the range 1/3 < v < 2/5 the composite 
fcrmion theory neglecting CF-LL mixing provides a very 
good description of the actual states obtained from exact 
diagonalization [7,8]. However, this approximation may 
become worse when the effective cyclotron energy is not 
too large; here the hybridization of CF states with higher 
CF Landau levels may become relevant. 

Given that the composite fermion system actually rep- 
resents a strongly correlated state of matter, the in- 
teraction energy for composite fermions in general in- 
volves two-, three-, and n-body terms. The most crit- 
ical assumption will be our neglect of three and higher 
body terms. Following Wojs and Quinn [18] we will as- 
sume that the interaction energy of the many composite 
fcrmion system can be modeled as a sum of pair-wise in- 
teraction. Furthermore, we will determine the pair inter- 
action from the system that contains only two composite 
fermions in the partially occupied CF Landau level (with 
the lower CF Landau levels fully occupied) . This assump- 
tion is valid only if the three and higher body terms do 
not cause a phase transition. 

There are indications that the higher order terms are 
small. The exact energy spectrum of the system with 
three composite fermions in the second CF-LL (where 
the only input is the Coulomb interaction between the 
electrons) has been compared with the energy spectrum 
that the pair-wise interaction predicts, and a reasonably 
good agreement is obtained [18]. The latter spectrum is 
obtained by diagonalizing the effective two-body interac- 
tion between composite fermions obtained by the method 
outlined earlier. However, no systematic study of the ef- 
fects of multi-particle interaction has been performed as 
a function of the number of particles. 

Another limitation of our work is its variational nature, 
which makes it as good as the trial wave functions for 
various states. For example, for FQHE states, we will 
use the standard wave functions for composite fermions 
[1]. These are extremely accurate in the lowest Landau 
level, but not so accurate for FQHE in higher Landau 
levels. Of course, it is also possible that the actual ground 
state has some completely new structure, not described 
physically by any of the trial states considered here. 

These caveats are meant to stress right at the outset 
that the results below are not "hard" predictions of the 
composite fcrmion theory, but are predicated upon sev- 
eral assumptions. This study is the first step toward a 
theoretical treatment of the interaction between compos- 
ite fermions, and it is hoped that future experimental and 
theoretical studies will shed light on the reliability of our 
present model. 



III. MODEL FOR INTERACTION BETWEEN 
COMPOSITE FERMIONS 

A. Pseudopotentials for CF-CF interaction 

It is known that a two-body interaction of fermions 
confined to some Landau level is fully characterized by 
its Haldane pseudopotentials [19], V mi which are the en- 
ergies of two particles in the relative angular momentum 
state m, IV'm): 



V m = 



(1>m\V\1> m ) 



(1) 
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Here Iq is the magnetic length, L m is a Laguerre polyno- 
mial and \ipm) is the two-body state of relative angular 
momentum to. For fully spin polarized fermions, only 
the pseudopotentials V m with odd integer values of to 
are relevant. In this paper, the Fourier transformation is 
defined symmetrically 
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where Jo is a Bessel function. In the equations below, the 
interelectron interaction, V(r) = e 2 /er, will be expressed 
in units of e 2 /el( )l and r in units of lg; the interelectron 
interaction then becomes V(r) = 1/r. The appropriate 
units will be put back at the end. 

Given the pscudoptentials we can construct the two 
body interaction in the lowest LL using the following 
formula. 



V(q) =2TrJ2V m L m (q 2 ). 



(3) 



Alternatively, we can directly write the Hamiltonian as 

ffW = EEw*-E E v m p%\i>). (4) 

i<j m i<j m—odd 

Here projects particles i and j onto a state of relative 
angular momentum to. This Hamiltonian is useful for 
exact diagonalization study. 

The pesudopotentials for composite fermions are de- 
fined similarly. For the pseudopotentials in the nth com- 
posite fermion Landau level, we consider the state in 
which the lower CF Landau levels are fully occupied and 
the (n + l) st CF Landau levels contains two composite 
fermions. V£ F is then defined as the energy of the state 
in which the two composite fermions are in the relative 
angular momentum to state. In other words, 
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where \4>m) 1S the wave function of composite fermions 
in the relative angular momentum m states. Note that 
the state of two composite fermions in the (n + l)st level 
is a many body state of electrons. For convenience, the 
energy is measured relative to E^ , the energy of the state 
in which the two composite fermions are far separated. 
(For finite systems the value of Eoo is chosen in the man- 
ner explained below.) 

We use the spherical geometry for our calculations. [19] 
Following the standard procedure for writing the wave 
functions for composite fermions [1] the wave function 
for the relevant state with angular momentum L is given 
by 



\1>Z F )=V L LL*l*t +L , 



(6) 



where $i is the wave function of one filled LL, is 
the wave function of the electron state in which n LLs are 
fully occupied and the (n+l)st LL contains two electrons 
in angular momentum L state: 



E 

mi ,m2 



< LM\lm 1 ;lm 2 > a+ +l rni a+ +l rn2 <5> n (7) 



Here $„ is the state with n-filled LLs, a+ +1 creates 
an electron in the l z — mi state in the (n + l)th LL, I 
is the angular momentum of a single composite fcrmion 
in the (n + l) st LL, and < LM\lmi; lm,2 > are Clebsch- 
Gordon coefficients. Because the energy depends only on 
the total angular momentum L, we work with M = 0, 
i.e., mi+m 2 = 0. The explicit form for the general wave 
functions of this kind is given in the literature [8]. Vlll 
is the lowest LL projection operator. 

In the spherical geometry, V m is defined as the energy 
of the state with angular momentum L = 21 — m [20]: 
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CF 
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This brings the definition of the pseudopotentials into 
conformity with the disk geometry, reproducing the disk 
pseudopotentials in the thermodynamic limit N — > oo, 
when the angular momentum of the lowest Landau level 
shell also diverges: I — > oo. To see this, it is helpful to 
note that in the disk geometry, the separation between 
the two particles increases with m. In the spherical ge- 
ometry, the smallest separation is obtained when both 
particles are put in the same l z state, which corresponds 
to the largest angular momentum L = 21 (not bother- 
ing about the Pauli principle here). Thus, the largest L 
corresponds to the smallest m. 

The integrals are performed using Monte Carlo; we 
evaluate V£ F for up to 50 particles. Of particular in- 
terest is the short range part of the interaction, which 



appears to converge to the thermodynamic limit fairly 
quickly with the number of particles, as expected. The 
long range part of the interaction is explicitly known, as 
discussed below. 

In order to minimize the computation time, we calcu- 
late several pseudopotentials within one Monte Carlo run 
by sampling with one wave function tp!~£(r), according 
to 
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CF 
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where (=) means that the two sides are equal upto an 
additive constant. Generally, the error increases as the 
angular momentum moves away from the sampling angu- 
lar momentum, and one must eventually use a different 
state for sampling; however, sampling with a single state 
works satisfactory for N < 40, where N is the total num- 
ber of particles. 

In Monte Carlo with composite fermions, all elements 
of the determinant are altered even if a single particle 
is moved, because of the strongly correlated nature of 
the problem. However, in a given step, the determinants 
that combine to give the state with a well defined L differ 
only in two columns. Rather than calculating each de- 
terminant entirely, we use an updating trick earlier used 
in Rcf. [21], which is a generalization of a method used 
for updating determinants that differ only in one column 
[22]. The method is explained in the Appendix. 

The pseudopotentials for the CF-CF interaction are 
shown in Fig. 1 of Ref. [17] and Fig (1). The explicit 
values are given in Tables I- VI. 



B. Inter-CF interaction in real space 

Now that the composite fermions of the fully occupied 
CF-LLs have been integrated out to produce an effective 
interaction between the composite fermions of the par- 
tially filled CF-LL, only the N composite fermions in the 
topmost partially filled LL will be explicitly considered 
in what follows. (Earlier, the symbol N was used for 
the total number of composite fermions, but below it will 
be reserved for the number of composite fermions in the 
partially filled CF-LL.) 

We thus have a system of fermions restricted to a cer- 
tain CF LL, with an interaction known through its pseu- 
dopotentials. It would be most natural to feed the pseu- 
dopotentials into an exact diagonalization routine and 
compute the ground state. However, such studies are 
not useful for some of the states that we wish to consider. 
The system sizes accessible to exact diagonalization are 
too small to capture the physics of compressible or stripe 
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states. (Also, the spherical geometry often used for ex- 
act diagonalization is not particularly convenient for the 
stripe phase.) 

We shall instead proceed by considering variational 
wave functions for several different kinds of states, and 
determining which has the lowest energy. As mentioned 
earlier, this approach has the shortcoming that we do 
not know if we have missed the true ground state, and 
even for a given state, we may not have a sufficiently ac- 
curate wave function. We will consider all of the states 
that have been investigated in the past in the lowest and 
higher Landau levels of electrons, and hope that the ac- 
tual ground state is one of them. Ultimately, only exper- 
iments can tell us for sure if that is the case. 

We will evaluate the energies of the some of the trial 
wave functions using the Monte Carlo method, for which 
we will need a real space form for the interaction be- 
tween composite fermions. Furthermore, it is convenient 
to write various trial wave functions in the lowest Landau 
level. So, we will map the problem of composite fermions 
at v* — n ± ^* into that of N fermions in the lowest LL 
at i>*. Such a mapping can, in principle, be carried out 
exactly, because the problem of fermions in any given LL 
with one interaction is equivalent to that of fermions in a 
different LL with another interaction, provided the Hal- 
dane pseudopotcntials for the two interactions are the 
same. Defining the problem in the lowest LL helps us 
to find the energies for various composite fermion states, 
because convenient techniques exist for writing the wave 
functions in explicit form in the lowest LL. 

For some of the states (paired state, Fermi sea) the 
calculation will be done in the spherical geometry, and 
the energy in the thermodynamic limit (N — > oo) will 
be estimated by an extrapolation of finite system results. 
For the stripe and bubble phases, the calculations will be 
done directly for the thermodynamic limit in the planar 
geometry. In this case, the spherical geometry is used 
only for the determination of the effective interaction be- 
tween composite fermions. 

In order to map the problem of composite fermions at 
v* = n±i>* into fermions at u* in the lowest Landau level, 
we need to find an effective interaction V e f / (r) which pro- 
duces the desired pseudopotentials in the lowest Landau 
level. A proper consideration of the long range part of 
the interaction is important for our purposes below, es- 
pecially in the context of stripes, the period of which 
is determined by a competition between the short range 
attraction and long range repulsion. Therefore, we treat 
the long distance part of the interaction explicitly in our 
calculations. Fortunately, the form of the interaction in 
this limit is known: it is given by the Coulomb inter- 
action between particles of fractional charge e/(2n + 1), 



i.e., 



(2n + l)" 



e 

do 



where r is the dimensionless distance between the com- 
posite fermions measured in units of Iq, the actual mag- 
netic length at v. At the effective magnetic field, i.e., 
the magnetic field corresponding to the filling factor v* , 
the natural unit for length is I*. We denote the distance 
measured in units of 1* by r*. The long range part then 
translates into: 

-5/2 



(2n+l)~ 2 l 
r* I* 



e 
el 



(2n+l)" 



e 
el 



(12) 



where we have used h/l* = (2n + I)" 1 / 2 . We model 
the full interaction between composite fermions at the 
effective magnetic field B* by the following form: 



V eff (r 



E 



*2i - 

Cir e 



+ 



el 



(13) 



(11) 



The power series, with properly chosen coefficients Cj, 
takes care of the short range part of the inter-CF in- 
teraction. The functions of the kind r* 2t e~ r are used 
because the pseudopotentials for such functions are ana- 
lytically calculable. Many other functional choices could 
work equally well; the mapping from pseudopotentials to 
a real space interaction in the lowest LL is one-to-many, 
and all real space interactions with the same pseudopo- 
tentials are completely equivalent in the absence of LL 
mixing. (Slightly different functions were chosen in Ref. 
[23]. There, the long rage Coulomb interaction was not 
considered explicitly, because that work dealt with the 
composite fermion Fermi sea, which corresponds to the 
limit n — > oo, so the last terms vanishes.) It is important 
to choose a real space form for the effective interaction 
that makes the interaction as smooth as possible; large 
oscillations of the interaction lead to larger errors in the 
Monte Carlo evaluation of the energy. 

The units for the effective potential defined above 
ought to be noted. If we replace e 2 /elo by e 2 /el*, then 
the interaction has the natural form for a problem at the 
effective magnetic field. We do our calculations with this 
interaction potential, and then obtain the final energies 
by a change of units. 

Next, we proceed to find values for the coefficients Ci 
so that the pseudopotentials of V e " reproduce the pseu- 
dopotentials V£ F calculated earlier from microscopic 
wave functions. Because we are considering fully polar- 
ized states, only the odd pseudopotcntials are relevant. 
In principle, it is possible to fix all of the Cj's to reproduce 
all V£ F, s exactly, but in practice, we will fix the first sev- 
eral (five to six) odd pseudopotentials only. That is suffi- 
cient for our purposes, given many other approximations 
that go into the calculation. The energy is determined 
predominantly by the short distance part of the interac- 
tion. Of course, at long distances, the above model inter- 
action automatically gives the correct behavior, thereby 
ensuring that it is quite accurate almost everywhere. 
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The energy of two composite fermion particles evalu- 
ated from microscopic wave functions gets contribution 
from two sources, their self energies and their interac- 
tion. The self energy can be identified by putting them 
far apart; because it is constant, it does not affect the 
calculation in any way. We find it convenient to subtract 
it out by shifting all V£ F by a constant so that the last 
pseudopotcntial (remember that there are only a finite 
number of pscudopotentials in a finite system) becomes 
equal to the pseudopotential of the long range part of the 
model interaction: 



-Boo = 



(ipr, 



| (2»+1) g 



IVWax) 



(14) 



where lipm^J) i s the two-body state with relative an- 
gular momentum m in the LLL. It is expected that for 
sufficiently large systems, the contribution of the gaus- 
sian terms to the last pseudopotcntial is negligible due 
to their rapid decrease with distance. 

We have matched the first five or six odd pseudopoten- 
tials by adjusting the same number of coefficients Cj. The 
values of Cj used below are quoted in Tables (VII) and 
(VIII). The comparison between the pseudopotentials 
V^ F , calculated from wave functions, and the pseudopo- 
tentials calculated from the real space model interaction 
are shown in Fig. 1 of Ref. [17] and 1, indicating that the 
real space interaction provides a good representation of 
the V^ F . The specific choice of which Cj are taken to be 
non-zero in our work (Tables VII and VIII) is dictated by 
the smoothness of the real space form of the interaction. 
We also note that the form of V e ^(r*) depends on de- 
tails, and can appear rather strange, with oscillations and 
even attractive regions, but proper results are guaranteed 
so long as it generates the correct pscudopotentials. 



IV. PRELIMINARY OBSERVATIONS 



composite fermions is trustworthy. (The caveat is im- 
portant; for example, three body interaction could play 
an important role, or when there are many composite 
fermions, the nature of the two-body interaction might 
change slightly from our model which considers only two 
composite fermion particles.) Nonetheless, the form of 
the interaction suggests that new FQHE states in the 
range 2/5 < v < 1/2 are either very weak or altogether 
absent. That is entirely consistent with experiments. 

For composite fermions in the lowest CF Landau level 
(u < 1/3), on the other hand, the interaction is strongly 
repulsive (Fig. 1, top panel), and new FQHE is ex- 
pected. That is also consistent with experiment. The 
2 CFs capture two additional flux quanta to become 4 CFs, 
which then fill Landau levels to produce FQHE at v = 
n/(4n± 1). There is presently experimental evidence for 
six members of these sequences [11]. These states can be 
interpreted either as the IQHE of 4 CFs or as the FQHE 
of 2 CFs; they are "pure," (as opposed to "mixed") in the 
sense that they contain only a single flavor of composite 
fermions. 

The CF pseudopotentials in Fig. 1 of Ref. [17] and 
Fig. (1) of this paper also indicate that the assump- 
tion of treating the interaction between many composite 
fermions as the sum of two body interactions is not un- 
reasonable. In a given CF-LL, Fig. 1 of Ref. [17] gives 
the interaction between two CF particles in an other- 
wise empty CF-LL, whereas Fig. 1 gives the interaction 
between two CF holes in an otherwise full CF-LL. The 
latter is a many body state of CF particles. However, the 
two interactions are rather similar both in shape as well 
as magnitude, which is what would be expected if only 
two body terms were relevant. 



V. VARIATIONAL STATES 
A. 4 CF FQHE / 4 CF Fermi sea 



A glance at the CF pseudopotentials reveals one of the 
most striking features of the inter-CF interaction: it is 
often attractive. The smallest energy is obtained at the 
smallest relative angular momentum, i.e., at the short- 
est distance. This is true for CF particles in the third 
or higher CF-Landau levels and for CF holes in the sec- 
ond and higher CF-LLs. This fact already appears to 
suggest an answer to one important question: Are there 
any new FQHE states between the fractions n/(2n +1)? 
New FQHE states here can appear if 2 CFs in a partially 
filled CF-LL capture additional vortices to form a FQHE 
state. However, capture of vortices requires a sufficiently 
strong short range repulsion. Given that the inter-CF 
interaction is actually attractive, any new FQHE in the 
range 2/5 < v < 1/2 appears unlikely (for fully polarized 
systems), to the extent that our model for interacting 



One possibility for composite fermions in the partially 
filled 2 CF-LL is to capture two additional vortices to turn 
into 4 CFs and then form an integral quantum Hall state 
with n' filled Landau levels, which corresponds to 



v* =n± 



filling for 2 CFs, and to 



2n'±l 



2v* ± 1 



(15) 



(16) 



for electrons. 

Taking 2 CFs as base particles, the wave function for 
their FQHE state is given by: 



2n' + l 



s 2 ^ 



(17) 
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where $„/ is the wave function of n fully filled LLs and 
Plll projects the state onto the lowest LL. Recall only 
the 2 CFs in the partially filled CF-LL are considered ex- 
plicitly here; the filled 2 CF Landau levels have been in- 
tegrated out in producing the effective interaction. The 
state at v* = n + 1/2 is obtained in the limit n' — > oo, 
where the filled Landau level state $„' becomes a Fermi 
sea. We expect these states to be energetically favorable 
when Vf !F is large compared to other pseudopotentials. 

The energies of these states can be calculated with the 
help of the real space interaction, following the method 
outlined in the literature. The spherical geometry will be 
used in our calculations, and the thermodynamic energy 
will be obtained from a linear, least squares fit of the 
finite system energies plotted as a function of 1/N. 



B. 4 CF pairing 

At half filled Landau level, pairing of composite 
fermions also becomes possible. A weak repulsion be- 
tween fermions is believed to lead to pairing of compos- 
ite fermions at half filling [12,13,16,23]. The second LL 
coulomb interaction between electrons falls within this 
class of interactions. The paired state is reasonably well 
described by a Pfaffian wave function [13], which, in par- 
ticular, has been shown to have lower energy than the CF 
Fermi sea at v = 5/2 [23]. Numerical studies on small 
systems have shown that this state has a fairly large over- 
lap with the exact ground state [16], although the over- 
laps are not decisive (in contract to the case for filled 
composite fermion Landau levels for which the overlaps 
are in excess of 99%). 

We ask whether or not the interaction between CFs at 
half filled CF LLs favors the paired state of 2 CFs. Here, 
2 CFs capture two vortices to turn into 4 CFs and then 
form pairs. With 2 CFs as base particles at filling factor 
equal to 1/2, we consider the Pfaffian wave function for 
the paired state: 



® Pf = $ 2 P/ [M] 



(18) 



where Pf [M] is the Pfaffian of the matrix M with 
components = (zt — z^y 1 defined as Pf[M] cx 

A [M 12 M 34 ...]. Pf [M] is a real space BCS wave function 
and so "Js can be viewed as a p-wave paired quantum 

2 

Hall state of 4 CFs. We calculate the energy of this state 
as discussed in the previous section. 



C. 2 CF stripes and bubble crystals 



stripes or bubble crystals. In this case, the base parti- 
cles, 2 CF, do not capture any additional vortices. We 
calculate the energies of these states closely following 
the Hartree-Fock formulation used earlier [24,25] for elec- 
trons in higher Landau levels; the only difference is the 
form of the interaction. We give here a brief outline of 
the calculation for completeness, which closely follows 
Rcf. [24]. 

Consider fermions in the plane of size L x x L y . In the 
lowest LL the eigenstates are 

Mx,y) = 1 e ^ e -^ +fc '°> 2 (19) 

where k £ jf- x (integer). Define the operator ^>(r) = 
Y]}. V'fe(r) afc where ak is the annihilation operator corre- 
sponding to the state \tpk >■ With the Fourier transform 
defined as 



and 



we have 



m = 



j- I d 2 q/(q)e^ r 



2tt 

L X Ly 



/(q)e* 



(20) 



(21) 



(22) 



The Fourier transform of the density operator, p(r) = 
^ (r)^f (r) , is given by 

= ^J d " r E^( r )^'( r ) e_<q ' r 4«fe' (23) 

k.k' 



With the following change of variables 
k + k' 



ly — k k, 



k 



k±=k T : f (24) 



the density operator becomes 



k 

1 ( K+fe+ ;2) 2 - 1 (x+k-llf-iq x x f 



--q l + ik q x l 



a' k+ a k _ 
a\ + a k _ 



(25) 



(26) 



For an attractive interaction, phase separation is a 
likely possibility, which, from our experience with higher 
Landau level physics, would lead to the formation of 



Define the interaction operator V: 



7 



V 



l - J d 2 vd 2 v' 



(2^) 3 
2L fr L, 



V(r,r>(r)p(r') 



^F(qM-q)p(q) 



= z^EE^)* 



exp 



q k ,k' 



4', o fc'_ a L°fc+ ( 2? ) 



In the Hartree-Fock (HF) approximation, the energy of 
the system is given by 



(V) 



L X Ly 

q,fco,fco 



E 



e 5<? 2 'o+ i ( fc o- fc o)gxio x 



(28) 



In terms of the operator A (q) = ^ fc e ~* fe9xi °(a|. a*;_), 
which represents the density of orbit centers, the energy 
is 



(V) 



(27T) 



2L T L, 



^[[/(g)-/^(^ 2 )]A(-q)A(q) 



q 



C/(|r. 



A(r)A(r') 



(29) 



where £/((?) = V(q)e~i g2l ° and U(r) is its Fourier trans- 
form. It is convenient to define Uhf{<i) = U (q) — IqU (qlg) 
where the first term corresponds to the direct interaction 
and the second to the exchange interaction. (V) in the 
above equation gives the CF-CF interaction energy, to 
which the CF-background and background-background 
interaction energies must be added to obtain a finite re- 
sult. Following Ref. ??, we define the cohesive energy as 
the energy difference measured from the uniform Hartree- 
Fock state 



E 



coh 



1 (27T) 3 

N 2L x L y X 

^[t/(<z)-? 2 C/(^)]A(-q)A(q) 

q#0 



(30) 



We can now calculate the cohesive energies of the stripe 
and bubble phases. In the stripe phase with the periodic 
length of stripe A, the orbit-center density is written as 



A( X , ) = J_ V 2sin(^) 



(31) 



In the bubble phase with the lattice constant for the 
hexagonal lattice, radius of the bubble R = \J : ^ L A b and 
area of a cell A = -^A 2 , the orbit-center density is 



A q ( o9 



(33) 



_ 2tt - _ 2ir 



where q = nej + me 2 , ei = ^y, e 2 - f^x - -ffi-y. 
This gives the cohesive energy 

^s^^^H'- (34) 

With our effective interaction the sums in stripe and 
bubble calculations converge quickly as we increase the 
cut off for q. Ref. [24], which deals with electrons in 
the higher electronic LLs, finds that the transition from 
the bubble to the stripe phase occurs at v n w 0.4, where 
v n is the electron filling in the (n + l) st Landau level. 
Somewhat surprisingly, we also find a transition at ^* w 
0.4 in all of our calculations. 

As far as the comparison of stripe and bubble phases 
is concerned, it is sufficient to know the cohesive energy; 
the energy of the reference, uniform HF state need not 
be evaluated. However, to compare with the Fermi sea, 
the FQHE state, or the paired state, we also need to 
know the energy of the uniform HF state, which must 
be subtracted from the full energy to obtain the cohesive 
energy. 

To calculate the energy of the uniform HF state we as- 
sume that the CF-CF, background-background, and CF- 
background interactions are all of the same form, charac- 
terized by V£ F . Since the direct interaction terms cancel 
each other the only remaining part is the exchange term 
of the CF-CF interaction. Thus the energy of the uniform 
state Eq is 



F U{0) x v 
hjQ — — x v 

For example in the case of the Coulomb interaction: 



(35) 



U(0) 



1 



E a = — x v = -- J dqez q X v = - , , - x u 

(36) 

which is a familiar result. For other interactions, we eval- 
uate Eq numerically. 



VI. RESULTS AND DISCUSSION 



where {q} = ?f x {• • • , -3, -2,-1, 1, 2, 3, • • •}. So the 
cohesive energy is 



E C oh 



2vl 



o 



Aq 



(32) 



The energies of various states are shown in Fig. 2 of 
Ref. [17], and Figs. (2) and (3) of the present manuscript. 
The lowest energy for the stripe or bubble states are de- 
termined by varying the period. We draw the following 
conclusions (subject to the validity of the model): 
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(i) FQHE: Within the lowest 2 CF Landau level, CF 
holes capture two more vortices to become 4 CFs. They 
show quantum Hall effect at is* = 1/3, which corresponds 
to 2 CF filling factor of v* = 1 — |, and electron filling 
factor of v = |. At i/* = 1/2, which corresponds to 
v = 1/4, the Fermi sea has the lowest energy. These 
results are consistent with the observations of FQHE at 
2/7 and CF Fermi sea at 1/4. 

(ii) Stripes: At v* = n + |, which correspond to 
v = (2n + l)/[2(2n + 2)], the stripe phase has the lowest 
energy (for n > 1). The stripe phase is obtained inde- 
pendently of whether we model = n+ i as 1/2 filling 
of CF particles on the background of n filled 2 CF Landau 
levels, or as 1/2 filling of CF holes on the background of 
n + 1 filled 2 CF Landau levels, which suggests that the 
result is robust. It should be noted, however, that the 
issues regarding the stability of the unidirectional charge 
density wave against a modulation along the length and 
also against quantum fluctuations, which have been in- 
vestigated in the context of electronic stripes in higher 
Landau levels [26], have not been considered here. 

The lowest energy for stripes is obtained for period 
A/l = 10, 28, and 34 for v = 3/8, 5/12, and 7/16. The 
period is rather large compared to that for the electron 
stripes in higher LLs (for which A/Iq is of order unity), 
which is not surprising because the interaction between 
composite fcrmions is rather weak, and also the difference 
between the densities of the FQHE states on either side 
is quite small. 

A transport anisotropy in higher electronic Landau lev- 
els, interpreted in terms of stripe formation, is observed 
at temperatures below ~ 50 mK [27]. The conditions for 
the CF stripes are more stringent. Estimates of the crit- 
ical temperature from the Hartree Fock theory are not 
quantitatively reliable, but noting that the effective in- 
teraction between composite fermions at v* = n+ 1/2 is 
roughly an order of magnitude smaller than for electrons 
at v = n + 1/2, as measured by the pseudopotcntials, 
we expect the critical temperature to also be similarly 
reduced. Also, the much larger period suggests the need 
for a high degree of density homogeneity. 

(hi) Bubbles: At v* = n ± | with n > 1, the bub- 
ble phase has the lowest energy. We remind the reader, 
however, that our trial wave functions for FQHE states 
that work very well in the lowest Landau level are not so 
good in higher Landau levels, which makes our calcula- 
tion somewhat biased against FQHE. 

To estimate at what filling a transition occurs from the 
bubble crystal to the stripe phase, we have determined 
their energies as a function of the filling factor, shown in 
Fig. (4) . The stripes are found to be stable approximately 
in the region 0.4 < v* < 0.6, outside of which bubbles 
take over. Overall, the phase diagram for various states 
of composite fermions, shown schematically in Fig. 5, is 
remarkably similar to that for electrons. For electrons: 
(i) FQHE occurs in the lowest Landau level; (ii) stripes 



are believed to be relevant in the vicinity of^ = n+ l/2 
for n > 2; (iii) bubble or Wigner crystal takes over in 
higher Landau levels for v = n + v' with v 1 < 0.4. The 
behavior for 2 CFs is quite analogous. 

VII. FINITE THICKNESS 

We have assumed until now that the electron layer 
width is zero. In actual experiments, the electron wave 
function has a finite extent in the transverse direction, 
which modifies the interaction pscudopotentials. The 
modified interaction has been obtained in a local-density 
approximation [28,29]. To estimate how finite thickness 
affects the results presented above, we have calculated 
the pseudopotentials for composite fermions in the sec- 
ond CF-LL and found (Fig. 6) that the value V% F , which 
is the largest pseudopotential, is reduced as we increase 
the density of the electrons. However, the change is not 
large enough to alter the previous results. As seen in 
Fig. (7), at v = 4/11, where the composite fermion fill- 
ing is v* =4/3, the bubble crystal phase continues to be 
most favorable among the ones studied. 

VIII. CONCLUSION 

In summary, we have considered theoretically the ques- 
tion of what new states of composite fermions are feasible 
as a result of the residual interaction between composite 
fcrmions. For this purpose, we have constructed a model 
for the CF-CF interaction, and studied various plausible 
states within a variational scheme. Our results suggest 
that fractional quantum Hall effect, stripe phase, as well 
as bubble crystal of composite fermions can all occur at 
various filling factors. 
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Appendix 

One of the most time consuming aspects of our Monte 
Carlo calculation is the evaluation of determinants. 
When two matrices have many common columns (rows), 
the determinant of one can be related to the other, which 
amounts to substantial saving of computation time in 
Monte Carlo. Here we describe a general method for it. 

Consider two n x n matrices A and A', which differ 
only in the (ii, Z2, • • • , i TO )th columns, where < m < 
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n. Our aim is to obtain the determinant of A' with 
minimum computation. The determinant of A' can be 
written as 

det[A'] = det[A] detfA^A'] (37) 

Since A' differs from A only at 12, ■ • • , z m )th columns, 
[A _1 A'] H is equal to 5 a when I <fc (ii, i 2 , ■ ■ ■ , i m ) 

Define Cj (j = 1, 2, • • • , m) to be n x n matrices whose 
components are same as those of A at jth column and 
zero elsewhere. Further, define an m x m matrix B such 
that [B] ab = [A- 1 Cj 3 ] 4a4b wher e a, b e {l,2,---,m}. 
Then the following is true. 

detlA^A'] = det[B] (38) 

This gives a simple relation between the determinants 
of A' and A: 

dct[A'] =dct[A]det[B] (39) 

Consider some examples: 
1. Updating of one column: 

dct[A'] = dct[A] dct[B] = detlAHA^C;];^ (40) 

2. Updating of two columns: 

det[A'] = det[A] det[B] 

= det[A]([A- 1 C' 1 ] nn [A- 1 C^] 4242 

-[A^CH^JA- 1 ^,,,) (41) 

3. Updating of three columns: 

det[A'] = det[A] det[B] 
= det[A] x 

([A 1 C 1 ]i 1 i 1 [A 1 C 2 ]i 2 i 2 [A 1 C 3 ]i 3 i3 

+ [A 1 C' 1 ] i2 i 1 [A 1 C2]i 3 i 2 [A 1 C' 3 ]i 1 i 3 

+ [A 1 C 1 ]i 3 i 1 [A 1 C 2 ]i 1 i 2 [A 1 C 3 ]i 2 i 3 

[A 1 C 1 ]j 2 j 1 [A 1 C 2 ]i 1 i 2 [A 1 C 3 ]j 3 j 3 

— [A C-J^JA C 2 ]i 2 i 2 [A Cgjjj^ 

-[A^C'^iA-^UAA- 1 ^]^) (42) 
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TABLE II. Pseudopotentials for CF-particles in the 2nd 
CF Landau level. 



m 




AE m 


1 


-11.755883 


0.0005593626 


3 


-11.741289 


0.0006092897 


5 


-11.756699 


0.0005888317 


7 


-11.754222 


0.0005821469 


9 


-11.754109 


0.0006913264 


11 


-11.755223 


0.0006993560 


13 


-11.755586 


0.0008570091 


15 


-11.755554 


0.0005545685 


17 


-11.756373 


0.0006021690 


19 


-11.756839 


0.0005796026 


21 


-11.756438 


0.0005914016 


23 


-11.756528 


0.0006283216 


25 


-11.756233 


0.0006583812 


27 


-11.756788 


0.0006599042 


TABLE III. 


Pseudopotentials for CF-holes 


in the 2nd CF 


Landau level. 






m 


Era 


AE m 


1 


-18.135544 


0.0009774497 


3 


-18.120577 


0.0009665822 


5 


-18.133069 


0.0009376358 


7 


-18.131464 


0.0008402889 


9 


-18.130383 


0.0009074839 


11 


-18.130536 


0.0009764340 


13 


-18.131314 


0.0009651683 


15 


-18.130754 


0.0010156039 


17 


-18.131685 


0.0009203913 


19 


-18.130155 


0.0009442089 


21 


-18.131228 


0.0010169967 



TABLE I. Pseudopotentials for CF-holes in the lowest CF 
Landau level. In this and the following tables, E m is the en- 
ergy of the full state, and AE m is the Monte Carlo statistical 
uncertainty. The energies are quoted in units of e /do- 





m 


E m 


AE m 


1 


-7.726983 


0.0005744120 


3 


-7.757470 


0.0007778533 


5 


-7.752651 


0.0010310604 


7 


-7.758716 


0.0010988030 


9 


-7.758750 


0.0012034221 


11 


-7.759422 


0.0012985310 


13 


-7.760764 


0.0012624134 


15 


-7.761201 


0.0012929297 


17 


-7.762082 


0.0016546843 


19 


-7.762933 


0.0016976433 
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TABLE IV. 


Pseudo-potentials for CF-particles in the 3rd 


CF Landau level. 




m 






1 


-18.942399 


0.0010950836 


3 


-18.932547 


c\ nm 01 e A o^n 
(J.UUlzlo4^7U 


5 


-18.927645 


0.0011055870 


7 


-18.935310 


0.0012656451 


9 


-18.934032 


0.001233024 


11 


-18.934417 


0.0012795043 


13 


-18.933864 


0.0013178322 


15 


-18.933060 


0.0012459279 


17 


-18.935017 


0.0013152929 


19 


-18.934135 


0.0013808266 


21 


-18.934547 


0.0015265486 


23 


-18.934336 


0.0012488644 



TABLE VIII. Parameters for our model interaction for two 
CF holes in an otherwise full CF-LL. The first, second and 
third CF-LLs are considered. All of the Ci's not shown here 
are set equal to zero. 



parameter 


1st CF-LL 


2nd CF-LL 


3rd CF-LL 


C'2 


2.08604 


-0.641674 


1.15620 


C6 


-0.0122001 


0.00240019 


-0.00699306 


ClO 


2.85410E-6 


-2.15587E-7 


1.38513E-6 


Cl4 


-7.63371E-11 


-1.27142E-12 


-2.49865E-11 


Cl8 


2.90620E-16 


1.48118E-17 


4.72082E-17 



TABLE V. Pseudopotentials for CF-holes in the 3rd CF 
Landau level. 



m 


E m 


AE m 


1 


-17.620082 


0.000983051 


3 


-17.612904 


0.000917553 


5 


-17.606845 


0.001063009 


7 


-17.614824 


0.001009564 


9 


-17.615345 


0.000978099 


11 


-17.614743 


0.000939712 


13 


-17.614587 


0.000958717 


15 


-17.612469 


0.000891311 



TABLE VI. Pseudopotentials for CF-particles in the 4th 
CF Landau level. 



m 


E m 


AE m 


1 


-16.691824 


0.00101485 


3 


-16.686846 


0.00105839 


5 


-16.682793 


0.00107136 


7 


-16.680260 


0.00108352 


9 


-16.687565 


0.00107323 


11 


-16.684489 


0.00111337 


13 


-16.685492 


0.00125616 


15 


-16.685188 


0.00095845 



TABLE VII. Parameters for our model interaction for two 
CF particles in an otherwise empty CF-LL. The second, third, 
and fourth CF-LLs are considered. All of the c;'s not shown 
here are set equal to zero. 



parameter 


2nd CF-LL 


3rd CF-LL 


4th CF-LL 


C2 


-1.57357 


1.75418 


-10.6790 


C6 


0.00817294 


-0.0114760 


0.0793665 


ClO 


-1.72040E-6 


2.62965E-6 


-2.49751E-5 


C14 


4.31056E-11 


-6.30827E-11 


9.87564E-10 


Cl8 


-1.80465E-16 


1.96255E-16 


-7.12056E-15 


C22 


0.0 


0.0 


9.29168E-21 
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FIG. 1. This figure gives the pseudopotentials for compos- 
ite fermion "holes" in the first, second, and third composite 
fermion levels. These are defined (up to an overall additive 
constant) as the energies of two composite fermion holes in the 
appropriate composite fermion Landau level, with the lower 
CF levels fully occupied. The systems are large enough that 
the pseudopotentials at the first few m (relative angular mo- 
mentum) values are well converged. Circles denote the pseu- 
dopotentials computed from the microscopic wave functions, 
with error bars showing the Monte Carlo uncertainty, and 
crosses are the pseudopotentials for the real space model in- 
teraction explained in the text. The pseudopotentials for CF 
"particles" in the second, third, and fourth composite fermion 
levels were given in Ref. [17]. 
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FIG. 2. The cohesive energy per particle for the CF Fermi 
sea (squares) and the CF paired state (circles) as a function 
of N, the number of CF "holes" in the topmost CF level. The 
thermodynamic energies are also shown for the CF stripe and 
bubble phases by a dash ( — ) and cross (x) on the y-axis, 
respectively. All energies are measured relative to the un- 
corrected uniform density state, explained in the text. The 
corresponding results for CF "particles" at v* = 1 + 1/2, 
v* — 2 + 1/2, and v* = 3 + 1/2, corresponding to v = 3/8, 
5/12, and 7/16, were given in Ref. [17]. 



13 



0.00 



-0.01 



-0.02 



.W 0.000 

>- 

a 

DC 
HI 



-0.005. 



LU 



0.005 



0.000 



-0.005 



v*=1-i 



v*=1+i 



----B-H3_ 



v*=2-4- 



v*=3-i 



v*=3+i 



0.005 



0.000 



-0.005 



0.005 



0.000 



-0.005 
0.005 



0.000 



-0.005 



0.00 0.02 0.04 0.06 

1/N 



0.00 0.02 0.04 0.06 0.08 

1/N 



-0.005 



-0.015 - 



bubble 

stripe 

v*=1-v* 




-0.005 



v*=2+v* 



v*=3+v* - -0.002 



-0.004 



-0.002 



-0.006 



-0.004 



0.000 



-0.008 



FIG. 3. The cohesive energy per particle for composite 
fermions at v* = n± 1/3, which correspond to electron filling 
factors v = 6 ^"f ± 2 » ^ or tne FQHE state as a function of N, 
the number of composite fermions in the partially filled com- 
posite fermion Landau level. This energy is shown by squares. 
The dash ( — ) and cross (x) on the y-axis show the cohesive 
energies of the stripe and bubble phases, respectively. 



FIG. 4. The cohesive energies of stripes (dashed line) and 
bubbles (solid line) as a function of the CF filling factor in 
various CF Landau levels. The striped state has lower energy 
in the vicinity of v* — n + 1/2. 
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FIG. 5. Phase diagram of various CF states, with bubbles 
(solid dots) , CF stripes (shaded region) and CF FQHE (solid 
line) as a function of the CF filling v* — n + u„. 



FIG. 6. Pseudopotentials for the CF-CF interaction in the 
2nd CF-LL for various densities for quantum wells of widths 
15 and 30 nm. The inter-electron interaction has been ob- 
tained in a local-density approximation. 
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FIG. 7. Same as in Fig. (3) for the modified interaction. 



1G 



